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Abstract. Enforcing exact conservation laws instead of average ones in 
statistical thermal models for relativistic heavy ion reactions gives raise to so 
called canonical effect, which can be used to explain some enhancement effects 
when going from elementary (e.g. pp) or small (pA) systems towards large AA 
systems. We review the recently developed method for computation of canonical 
statistical thermodynamics, and give an insight when this is needed in analysis of 
experimental data. 
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1. Introduction 

Statistical thermal models are widely used in the analysis of relativistic heavy ion 
collisions (see e.g. |Q and references therein). In these models, the final stage 

of inelastic interactions between hadrons, the chemical freeze-out, is characterized by 
the parameters of stationary maximum entropy ensemble. 

For large system, it is appropriate to choose the most straightforwardly applicable 
ensemble, the grand canonical (GC) one, which is parametrized by temperature T 
and chemical potential Hi for each averagely conserved charge i. The GC ensemble is 
defined in the large volume limit, so the volume parameter V is an extensive coefficient 
in the expressions for such quantities as mean particle numbers. For small systems, 
such as pp or pA reactions, the charge conservation laws must be handled exactly, so 
the ensemble to be chosen is the canonical (C) one. Then, the V independent fugacities 
due to average conservation laws do not appear, but are replaced by nonlinearly 
V dependent canonical chemical factors, which can account for many enhancement 
patterns going from pp to pA to AA systems B 110- 

In this paper, we quote the newly developed computational method for relativistic 
B, S, Q canonical hadron ensemble, which enables us to reach previously unattainable 
canonical model results, up to baryon number 0(100) whereas with old methods it 
was only possible to reach B ~ 20 |^] with very long computation times. Once 
this is achieved, we can give an answer to question whether the full B, S, Q canonical 
modelling is needed, or is the S canonical {B, Q grand canonical) or GG approximation 
accurate enough. 


2. Methods 


In order to pass from the usual grand-canonical partition function Zqc to the one 
fulfilling exact internal symmetries, we employ a well-known group theoretical method, 
first introduced by Gerulus i I, 0- By denoting the set of conserved quantum 
numbers by {Ci}, the canonical partition function can be obtained by using a 

projecting operator onto the conserved quantum numbers. In the case of N internal 
symmetries of type U{1), the projection takes the form: 


Z{cMT.V) 



Zgc{T, V, {AcJ). 


( 1 ) 


where (fi S [0, 27r) is a U{\) group parameter and a Wick-rotated fugacity factor 
Aci = is introduced, for every charge Ci. 

In heavy ion reactions, the relevant set of conserved charges is {Ci} = B,S,Q, 
namely baryon number, strangeness and electric charge. Applying this set to the 
previous equation yields a triple integral, which is computationally extremely time 
consuming for B > 5. In p]| we have eliminated analytically the baryon integration, 
and have obtained a very efficient way to calculate the full canonical {B, S, Q) hadron 
thermodynamics. The applied form of partition function reads: 

ry n‘2'K p27T 

Zb,s,q{T,V) = ° / d(j)s j d(/)Q 
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( 2 ) 
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exp 


2 X! ^OsiSi(j)s + Q^(t)Q) 
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where Zq denotes the partition function for hadrons carrying no relevant charge and 
/ is the modified Bessel function. The sum is over mesons, and zj is the one- 
particle partition function. The vj above is defined as zf where 

the summation runs over baryons but not antibaryons. Now, we are left with a double 
integration only which can be then performed numerically with no major problem. 

The mean particle numbers of primary hadrons (i.e. those directly emitted from 
the hadronizing source) (Ni) are obtained from the equation (||) |^, by taking the 
derivative of the canonical partition function with respect to a hctitious fugacity : 


(iV,) = A, 


d\nZB,s,Q{T, V) 


\i = l 

ZB-Bi,S-Si,Q-Qi{T,V) 1 
-Z,- , 


Zb,s,q{T,V) 

where the quantity ZB-Bi,s-Si,Q-QilZB,s,Q is called chemical factor 1^. In the 
large volume (thermodynamical) limit, this expression becomes the GC one, i.e. 
(iV,) = Ai-AfU§‘zi. 


3. Numerical results 

In this paper, we perform no actual htting procedure to the experimental data, but 
only use some suggestive thermal parameter values from our previous work ||^ . This is 
because our main motivation is to show systematically the canonical effects in different 
systems in order to give the reader an insight on applicability of different thermal 
approaches, namely full canonical, strangeness canonical and grand canonical. 

In the following, we fix the temperature and baryon density in several cases. 
Strangeness and charge to baryons ratio are fixed according to initial nuclear 
composition. These conditions lead to fixed fugacities and particle ratios in GC 
formalism. GC approximations are then to be compared with full canonical and 
S'-canonical counterparts, which are nonlinearly dependent on net baryon content. 

In figure we show the canonical kaon enhancement in the conditions relevant to 
GSI energies, 1.7 GeV per nucleon in SIS AuAu collisions. The theoretical production 
ratio K"*" over number of participating nucleons is shown to point out the canonical 
strangeness effect. Note that we have been able to compute the full chemical factor 
even for B = 400. The ratio increases very slowly towards the GC limit, thus 
strangeness must always be handled canonically when analysing SIS results. The 
relative difference between full canonical and strangeness-canonical results is 13% for 
B = 2, decreasing to 4% for B = 10 and to 1% for B = 40. It must be pointed out 
that the above results are calculated using an isospin symmetric initial configuration 
whereas Z/A ~ 0.4 for gold nucleus. However, this variation essentially gives no change 
on the relative differences above, although the ratio decreases naturally due 

to the initial neutron excess. 

In Au-Au collisions at AGS energies (11.6 A GeV momentum), the chemical 
freeze-out temperature is found to be around T = 120 MeV and the baryon density 
close to the normal nuclear density 0 . In hgure § we plot the calculated K '^with 
results obtained in peripheral to central Au-Au, Si-Au and Si-Al collisions. It can 
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be seen that all the way up from i? = 2 to i? = 60 the full canonical and strangeness- 
canonical results are essentially the same. All central reaction results lie on the region 
where the canonical effects are negligible and the GC formalism applies. Strangeness 
enhancement in Au-Au system does not definitely look like being of canonical origin, 
whereas Si-Au and Si-Al follow roughly the curve at the normal nuclear density (not 
shown in the figure). The above argument also applies for the K“/iVpart enhancement 
pattern Q. 

The best fit temperature for the multiplicities measured by NA49 experiment 
at SPS in central Pb-Pb reactions at 158 A GeV is found to be T 160 MeV Q. 
In figure we show the strange baryon enhancement by using a baryon density 0.3 
fm“^. Although the baryon density in our analysis [|] was found to be ub ~ 0.2 
fm“^ this larger value is used in order to probe the largest reasonable canonical effect. 
All results in figure ^ are normalized to the baryon multiplicities per participant at 
the point = 2. This choice reveals the slight difference between the results 

obtained using the S'-canonical approximation and the B, S, Q-canonical calculation. 
An interesting feature here is the fact, that the S'-canonical approximation leads to 
an overestimation of the canonical effect. 
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Figure 1. Canonical enhancement of kaons as a function of number of 
participants in the conditions relevant to GSI energies. 
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Figure 2. Theoretical K'^/A^part curves at fixed temperature to two different 
baryon densities shown along with AGS experimental results ||ig|]. Curves are 
strangeness canonical while squares are full canonical results 



Figure 3. Strange baryon enhancement in the conditions relevant to the SPS 
energies. Hadron multiplicities are normalized to results at H = 2. Curves are 
strangeness canonical while crosses are full canonical results. 














